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\^ ", We consider the infinite-range deterministic spin models with Hamil- 

tonian H = E£=i Ji 

jCFiCfj, where J is the quantization of a chaotic map 
£>v of the torus. The mean field (TAP) equations are derived by summing 



in 
o 

On 



the high temperature expansion. They predict a glassy phase transition 
at the critical temperature T ~ 0.8. 



PACS N. 64.60 Cn, 64.60 Fr, 64.70 Pf 



A class of infinite-range deterministic Ising spin models with glassy behaviour in nu- 
merical simulation has been recently identified flU, f2j, ||: unlike the random coupling 
case, however, the mean field equations (hereafter the PP equations) are different 
from the standard TAP ones || and their linearization does not determine a critical 



temperature for the glassy transition. 



a 

o . 

The most interesting representatives of the class are the sine (equivalently, cosine) 



models, with Hamiltonian 

5-H N 

# = -1/2 J2 JijWi (!) 

where the coupling matrix J is (twice the uppermost left block) of the discrete sine 
(cosine) Fourier transform 



2irij 

V2N + 1 "~ \2N+1 



In fact, the ground state of H can be explicitly computed || if 2A^ + 1 is prime with 
N odd: this leads to detect a first-order crystalline phase transition (at a higher critical 
temperature than the glassy one fHH), whose relevance on the glassy behaviour of the 
system is discussed in §3. In turn, the matrix J coincides |J with (the imaginary 
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part of) the unitary propagator quantizing the discrete dynamics generated by the unit 

fO -1\ 

symplectic matrix S — I acting as a Hamiltonian map over the 2— torus T 2 (the 

U / 

operator quantizing a Hamiltonian map of the torus is a N x N unitary matrix |7]||, 

N being the inverse of the Planck constant :Q in this context the thermodynamic limit 

N — > oo is formally equivalent to the classical limit). This algebraic identity suggests 

as more natural candidates for detecting deterministic glassy behaviour the coupling 

matrices defined by the quantization of hyperbolic maps over T 2 . The corresponding 

discrete dynamical systems are indeed chaotic, while S generates a period-4 one. Here 

fa b\ 

we consider the matrices A = with a = d = 2g,b = l,c = 4g 2 — l,g G N, 

Vc d) 

which admit two positive eigenvalues Ai > 1, A 2 < 1. The quantization of A is ]7|-|j| 
the unitary N x N matrix 

V(A) jk = CnN- 1 ' 2 exp [(2ttz/A0 (gf ~jk + gk 2 )} (3) 

with | Cat | = 1. The present models are thus defined by the Hamiltonians Ha{(?) = 
^2J(A) jk cr j a k with J (A) = Re[V(A)), i.e.: 

J(A) 3k = C N N- 1 ' 2 cos [{2m/N){gj 2 - jk + gk 2 )} (4) 

Now V(A) = V T (A) = V~\A) whence a(Re[V(A)\) = Re[a(V(A))\. Hence the 
equistribution [10| of the spectrum a(V(A)) over the unit circle for iV — > oo with an 



eigenvalue at 1 for all N implies the same properties for a (J (A)) in [—1, 1]. 

Our claim is that these deterministic models behave, at least in mean field approx- 



imation, more closely to the random ones, such as the SK one We will indeed 



The physical intuition is: the phase space has volume 1, and can accommodate at most N 
quantum states of volume h, so that Nh = 1. 
2 This holds only if a, b, c, d are as above, and motivates our choice of hyperbolic matrices. 



Deterministic spin models with a glassy phase transition 3 

compute in closed form (at the thermodynamic limit iV — > oo) the Gibbs (i.e., magne- 
tization dependent) free energy The result is: 

N 

(3® = 1/2 {(1 + mi)]n(l + ™<)/2 + (1 - ^i)ln(l - m,)/2} - 
- /3/2 J2 JMijimmj - NG{(3{\ -q))- G{(3) = (3 2 /(8 + A(3 2 ) (5) 

n 

where raj : Z = 1, . . . , N is the magnetization at site I and q = 1/N ^ m 2 the Edwards- 

z=i 

Anderson order parameter (as in [J3J, does not depend on the particular choice of 



J(A)). The stationarity condition of the Gibbs free energy |12| yields the mean field 



equations of the model (believed exact because of its infinite range) 

Qi = tanlT 1 rm + 2pG'(p(l - q))rrn - (5 ^ J(A) ijmj = (6) 

j 

Unlike the PP ones (but like the TAP ||) these equations can be solved by lineariza- 
tion near q = 0: any eigenvector for the eigenvalue 1 of J (A) yields a solution if 
2/? 2 /( 2 + I3 2 ) 2 -0 + 1 = 0. Its only positive zero (3 ~ 1.25 determines the critical tem- 

N 

perature T ~ 0.8. The phase transition is glassy because, as we will see, 1/iV^ raj — > 

i=i 

as N —>■ oo. This fact and ([5]) implies that all such solutions have the same free en- 
ergy. The thermodynamics is thus independent of their number, the multiplicity of the 



eigenvalue 1, which depends in a sensitive way on the integer iV \W . 

To get the mean field equations (^) we resum the high temperature expansion by 
the same procedure of M. The function G(/3) has a simpler form in this case because 
there is only one class of non vanishing diagrams at the thermodynamic limit. 
Consider the Helmholtz free energy F{(3): if Z(/3) = ex P PHa((t)/2 = exp —/3F(/3) 
is the partition function (at site-dependent magnetic field hi = 0) we have ]13] 

e-P F = (2/vr) Ar / 2 /det^(/?J) / exp{((2/3J)-V, <P) + $>gcosh(& + h)} ^0k=o = 

Jr n - 

2 n tt- n / 2 / exp{- < x,x > +^ log cosh [^£((J 1/2 ka;J } d N x (7) 

Jrn I . =i LV - jj 
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The high temperature expansion for /3F(/3) is generated out of the integration of the 
expansion of exp log cosh x . The well known diagrammatic representation of its n-th 



order term is obtained [13] by drawing all diagrams with n links, 2<j + l<n + l 
vertices and no external legs, whose individual contribution is 

1. For any link between two consecutive vertices I ^ k a factor (3Jik', 

2. For any vertex with m links the cumulant u m , i.e. the m— th coefficient of the Taylor 
expansion of log cosh x ( (T^j, p. 414); 

3. Any diagram has to be divided by its order of symmetry. 

The contribution of each individual diagram D at order (5 n is indeed 

\D\ = U(D)S(D)- 1 £ J n,r 2 -J%r 3 ---J% w j = l,...,n (8) 

n,...,r j+1 

Here j + 1 is the number of vertices, «i + . . . + otj — n; n — j + 1 = n — 1, . . . , 1 is the 
number of loops, > 1 the number of links between consecutive vertices, S(D) the 
symmetry factor, and U(D) = iti (011)1*2(0:1 + a^) " ' ' u j( a j-i + a j) u j+i( a j)- Remark 
that r n+ i = 7*1 for j = n. We now verify that, at the thermodynamic limit iV — > 00: 

(a) All diagrams for n = 2p + 1 vanish; 

(b) For n = 2p the only surviving diagram D p is the one with p + 1 vertices, p loops 
and two links between consecutive loops: 

00""O : w = i-iy- 1 ^- 2 



To prove (a) and (b), first recall the basic estimate fulfilled by the Gauss sums [14 



N 

]T exp (2m/N)g(s h . . . , < CN 1 ' 2 (9) 

Sl,...,Sl=l 

where g is any quadratic form in the / integers si, . . . , sj with integer coefficients and 
C a constant independent of g and N. By (|3|) and (§), the sum in @ amounts for 
any fixed 1 < j < n to 2 n Gaussian sums over r±, . . . ,Vj divided by N n l 2 . We now 
argue that all these Gaussian sums are, once divided by N n ^ 2 , of order 1 or less except 
one: that with n = 2p, oi = . . . = a p — 2, and summand independent of 7*1, ... , r p +i, 
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generated by the constant term in the expansion of J 2 k Th 1 ■ This sum is clearly 

k 

equal to jV p+1 jV~ p 2" p = N2~ p . Any other (divided) sum is down by at least N 1 : each 
summation over q indices is estimated by N q l 2 if the summand is a Gaussian and by N q 
if the summand is constant, which case shows up whenever a, is even for at least one 
i; on the other hand any power of aj in excess of one increases by one the number of 
loops and thus reduces by one the number of vertices and hence of summation indices. 
Thus for n = 2p + 1 the most divergent sums behave as N p+1 ^ 2 : this happens either for 
oij = 1 for all i so that j = 2p, or for ol\ = . . . = ct v -\ = 2, a p = 3 (and permutations 
thereof) so that j — p + 1. The prefactor jV _(i,+1 / 2 ) yields the estimate 0(1). For 
n = 2p the second most divergent sums have the same behaviour: again they take 
place for aj = 1 for all % and j = 2p — 1 or for ot\ = . . . = a p = 2 with the summand 
depending but on two indices (up to permutations). Now the only surviving D p has 
two vertices with two links (the extrema) and p — 1 vertices with four links (all the 
remaining ones). Since 112 = 1,W4 = —2, U(D P ) = (— 2) p . Moreover the symmetry 
factor S(D P ) is 2 • 2 P (to account for the interchange of the external vertices and of any 
pair of links between the internal vertices). Hence U(D P )/ S(D P ) = (— l) p_1 /4 and by 
(^D we have \D p \ = (— l) p_1 • 2~ p ~ 2 . Summing up we get the Helmholtz free energy: 

00 

- PF(0) - Mn2 = N ^(-l)^ 1 ■ 2- p ~ 2 p 2p = Np 2 /(8 + A(3 2 ) = NG{(3) (10) 

P =i 

As in [|j], expanding —f3F(f3) to first order we recover the SK Helmholtz free energy 
Mn2 + (3 2 /AJ2ij Jfj because here £ij J 2 d = N/2 + 0(1) for N large. 

To obtain the Gibbs free energy we have to perform the Legendre transform 

&(P,mi) =max hi [F(P,hi) -J2 h i m i\ ( n ) 

i 

out of the (hi dependent) expansion of (0) in powers of /3. To this end we simply 
take over the Parisi-Potters argument because we are summing over a subclass of 
the " cactus" diagrams considered in [[| within the same assumption of self-averaging, 
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N 

namely m 2 k = iV _1 E m ? + o(l) as iV — ► oo. This yields (|5|), whence the mean field 
1=1 

equations (|6]) and the critical temperature. 

N-l 

It remains to prove that N^ 1 E m l k — > as iV — > oo if J(A)m k = m k . Any such 
eigenvector defines indeed a critical (staggered) magnetization distribution; further- 
more the vanishing of the total magnetization represents the necessary condition both 
for the glassy nature of the transition as well as for self-averaging. Setting O = 

N-l 

N-V 2 (l, 1 . . . , 1) we have N' 1 E m[ = (m k ,4> ) = (m k ,V(A) l 4> ) VI E N. On the 

1=0 

other hand [JTOj V(A) tends (weakly) to a unitary operator on L 2 (—n, ir) with Lebesgue 



spectrum on the unit circle as iV — ► oo. This entails [TH] lim lim (m^., V(A) O ) 

£-^oo iV— >oo 



lim (m fc , O ) = 0. 

iV^OO 

Moreover by (§) the specific free energy / = is the same for all eigenvectors m k . 

the third term depends only on q, the second is —(5q/2 because Jm k = m k and 

N 



2iV -1 E {(! + mi)ln(l + m<)/2 + (1 - mi)]n(l - mi)/2} = g + 0(g 3/2 ) 

i=l 



since g is by definition small near the critical point. We conclude with two remarks, 
(i) The staggered magnetizations can be explicitly computed for some particular values 
of N. There is indeed |L0| p(N) G N such that J(A) P = I d , and for "most" sequences 
of values of N one has p(N)/N = M, M < oo. Under these conditions the eigenvalue 
1 of J (A) has multiplicity M + 1. A first corresponding (norm e) eigenvector is |J 

2 

mi = cos (2iri/N)kl 2 /y/eN where k = 3 (modiV), while 

p-i 



r 

m, 



1 /VtP E c « cos (2iri/N) (a r f + b r s l)/VN, r = 2, . . . , M + 1 



s=0 



where |c£| < 1 and < a r s ,b r s < N — 1 are integers ( ||, Formula (4.25)). By (H) one 

AT-1 

checks directly that iV -1 E m [ ~~ * f° r a U eigenvectors. 

1=0 

(ii) The magnetization just below the critical point has the square-root behaviour of 
a second order transition. Compute indeed the second order expansion 5 2 Qi near j3 c 
and m* as in |5J , putting (3 = (3 C + A/5 and rrii = m* + Smi where m* is any staggered 
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magnetization vector and (5m, m*) = 0. Neglecting the term of order A/3 2 , taking 
the scalar product with m* and dividing by Nq the equation (m*,5 2 Q) = we get 
a\\5m\\ 2 + 7 A/3 = with a = -2(3*G" (J3 C ) > 0, 7 = -1 + 2G'((3 C ) + 2(3 c G n ((3 c ) < 0, 
whence the assertion. 
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